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Aspects  of  Cardioid  Processing. 
D.T.  Hughes 


Executive  Summary:  Traditional  active  towed  array  sonar  systems  suffer 

from  the  problem  of  ambiguous  returns:  in  which  a  sonar  operator  cannot  ascer¬ 
tain  whether  a  possible  target  is  placed  to  the  port  or  to  starboard.  Although 
operational  methods  exist  for  overcoming  this  problem  they  are  time-consuming 
and  cannot  usually  be  used  for  a  single  ping. 

Several  NATO  nations  are  investigating  ‘triplet  technology’  to  allow  port- 
starboard  (left-right)  discrimination.  SACLANT  Undersea  Research  Centre 
has  purchased  a  broad-band  triplet  array  to  investigate  the  use  and  implemen¬ 
tation  of  this  important  and  relevant  technology. 

This  report  describes,  in  detail,  two  modes  of  operation  of  the  triplet  array. 
These  methods  are  investigated  for  SNR  performance  and  compared  for  corre¬ 
lated  and  uncorrelated  noise.  The  latter  is  expected  to  play  a  crucial  role  in 
this  technology  for  which  intra- triplet  hydrophones  are  closely  spaced.  Both 
methods  are  shown  to  give  excellent  results  on  real  data,  with  a  typical  Left- 
Right  suppression  of  greater  than  15dB,  allowing  unambiguous  determination 
of  target  bearing  as  well  as  reverberation  suppression. 
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Aspects  of  Cardioid  Processing. 
D.T.  Hughes 


Abstract: 

It  is  well  known  that  conventional  active  towed  arrays  suffer  from  ‘The  Am¬ 
biguity  Problem’  in  which  it  is  impossible  to  distinguish  returns  from  port  or 
starboard.  Although  operational  methods  for  overcoming  this  problem  exist 
they  are  time-consuming  and  cannot  be  used  for  a  single  ping.  One  proposed 
method  of  coping  with  these  difficulties  is  to  use  a  triplet  array  in  which  the 
direction  of  arrival  of  a  signal  can  be  ascertained  on  a  single  ping. 

This  report  details  a  broadband  algorithm  which  has  been  developed  and  is  cur¬ 
rently  being  used  at  the  Saclant  Undersea  Research  Centre.  The  algorithm  was 
specifically  designed  to  be  used  on  a  delay-type  beamformer  in  which  shading 
cannot  be  applied  at  the  beamlevel.  A  second  algorithm  employing  the  full- 
flexability  of  a  3d  beamformer  is  also  described  and  the  two  algorithms  are 
compared  for  their  performance  in  terms  of  left-right  suppression  and  SNR  for 
which  analytical  expressions  are  computed.  It  is  shown  that  both  of  the  al¬ 
gorithms  have  desirable  but  different  features.  The  algorithms  are  applied  to 
data  taken  from  the  trial  BACCHUS ’98  and  the  results  obtained  are  used  to 
help  validate  the  theoretical  results  and  to  give  an  indication  of  the  perfor¬ 
mance  which  is  achievable  for  the  triplet  system.  We  show  that  for  the  Centre 
algorithm  a  typical  achievable  left-right  suppression  is  greater  than  15dB. 
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Introduction 


Several  NATO  nations  have  an  active  research  and  development  program  investi¬ 
gating  aspects  of  Left-Right  suppression/discrimination  (including  UK,  France,  The 
Netherlands  and  Germany).  [1,  2,  3] 

The  present  operational  means  of  discriminating  whether  a  particular  sonar  return 
is  from  port  or  starboard,  in  effect  necessitates,  tracking  a  contact  over  many  pings 
from  a  manoeuvering  platform.  The  techniques  which  are  currently  under  investi¬ 
gation,  however,  discriminate  on  a  single  ping. 

Various  means  exist  for  carrying  out  this  discrimination  including: 

i)  directional  transmitters;  in  which  different  pulses  are  transmitted  to  port  and 
starboard. 

ii)  Twin  (or  multiple)  arrays;  in  which  two  or  more  line  arrays  are  towed  parallel 
with,  in  the  twin  case,  a  spacing  of  approximately  quarter  of  a  wavelength  of  the 
highest  frequency. 

iii)  Cardioid  technology.  This  can  be  further  broken  into  two  areas  of  interest.  It  is 
possible  to  produce  directional  transducers  allowing  an  array  to  be  made  which  is 
actually  formed  from  a  left  and  a  right  looking  array  ‘hard-wired  in’.  The  alternative 
approach  is  to  produce  an  array  made  up  of  triplets  of  sensors  with  all  channels 
digitized  and  sent  back  to  a  central  processor  so  that  beamforming  and  hydrophone 
sensitivity  can  be  controlled  in  real  time. 

It  is  this  last  technique  which  we  shall  investigate  in  this  report.  Figure  1  shows 
an  artists  impression  of  such  an  array.  Each  triplet  (for  example  that  marked  ABC 
on  the  figure)  provides  3  channels  of  acoustic  data.  The  individual  triplets  may 
be  rotated  away  from  the  vertical  but  the  degree  of  rotation  will  be  known  fairly 
accurately  due  to  roll  sensors  positioned  along  the  array.  (We  discuss  in  annex  E 
the  effect  of  errors  in  this  measurement).  Generally  the  array  is  constructed  in  such 
a  way  as  to  stop  the  array  from  rotating  too  severely  whilst  being  towed. 

S  ACL  ANT  Research  Centre  obtained  a  broadband  cardioid  array  in  1999.  This 
report  details  the  proposed  algorithm  for  use  in  the  Centre’s  system.  It  goes  on  to 
discuss  alternative  algorithms,  and  details  certain  important  operational  features  of 
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a  cardioid  array. 

The  philosophy  of  this  report  is  to  give  the  main  algorithm  and  indication  of  results 
in  the  body  of  the  text  whilst  leaving  the  detailed  mathematical  investigation  to  the 
annexes. 

In  section  2  we  detail  the  algorithm  which  is  to  be  used  in  the  Centre’s  array.  This 
algorithm  needs  to  be  broadband  and  implement  able  on  a  delay  beamformer.  We 
refer  to  this  algorithm  here  and  throughout  the  report  as  the  ‘Centre-algorithm’. 
The  general  approach  is  to  do  the  beamforming  in  two  phases.  The  first  is  done  at 
the  triplet  level  using  signal  processing  to  form  a  cardioid  responses  before  doing 
the  processing  in  the  azimuthal  plane. 

In  sections  3  and  4  we  briefly  discuss  the  philosophy  and  operating  parameters  of 
the  BACCHUS’98  trial  as  well  as  giving  some  typical  results  using  the  algorithm  of 
section  2.  In  section  5  we  draw  conclusions  and  make  recommendations  based  on 
the  main  body  of  the  text  and  from  the  annexes. 

In  annex  A  we  show  how  the  algorithm  of  section  2  can  be  obtained  as  an  approxima¬ 
tion  of  an  optimum  (in  a  least  mean  squares  sense)  beamformer  with  the  constraint 
of  real  weights. 

In  annex  B  we  discuss  one  possible  alternative  algorithm  which  may  be  employed  on 
a  triplet  array.  We  refer  to  this  as  the  ‘3d-algorithm’  since  it  treats  the  array  as  a 
whole  rather  than  breaking  it  into  the  individual  triplets  as  for  the  Centre- algorithm. 

In  annexes  C  and  D  we  compare  the  Centre-algorithm  and  the  3d-algorithm  in  terms 
of  SNR  and  left-right  suppression  statistics. 

In  annex  E  we  briefly  consider  how  the  algorithms  are  affected  by  system  errors  and 
approximations . 
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Figure  1  Artist’s  impression  of  Cardioid  Array .  The  triangle  ABC  represents  a 
triplet  of  the  array  ivhich  provides  three  channels  of  acoustic  data: 
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Broadband  Cardioid  Processing 


The  array  obtained  by  SACLANT  Undersea  Research  Centre,  as  well  as  being  di¬ 
rectional  has  been  designed  to  be  used  for  broadband  research,  with  a  bandwidth 
of  approximately  2.5kHz.  Consequently  any  proposed  algorithm  must  be  amenable 
to  broadband  processing.  Also  the  algorithm  should  be  implementable  in  a  time- 
domain  delay,  beamformer  rather  than  in  the  frequency  domain.  Furthermore,  due 
to  the  nature  of  the  Centre’s  real-time  system,  when  this  research  was  being  car¬ 
ried  out,  shading  cannot  be  applied  at  the  beam  level.  Subsequent  upgrades  to  the 
system  have  removed  this  restriction  but  the  algorithm  detailed  here  will  be  con¬ 
tinued  to  be  used  in  the  forseeable  future.  Such  problems  have  recently  received 
some  theoretical  consideration  in  the  literature  [5].  In  annex  B  we  discuss  another 
algorithm  which  although  in  so  me  ways  is  superior  to  the  Centre  algorithm  requires 
a  frequency  domain  implementation. 

It  is  also  worth  mentioning  at  this  point  that  other  algorithms  can  be  envisaged  and 
the  exact  form  that  we  expect  an  algorithm  to  take  is  dependent  on  what  is  desired 
from  the  final  system.  For  instance  in  the  two  algorithms  investigated  in  this  report 
we  do  not  use  the  ‘vertical’  extension  of  the  triplet  which  would  allow,  for  instance, 
a  null  to  be  driven  towards  the  surface  to  allow  suppression  of  surface  noise.  Such 
research  although  interesting  and  worthwhile  is  beyond  the  scope  of  this  report. 

In  what  follows  we  assume  that  there  are  Nc  triplets  (3JVC  hydrophones).  The 
hydrophones  are  numbered 


first  triplet  last  triplet 

(^11,  ®21,«31, . X1 N,  X2N i  x3n) 

(yiUV2U  2/31, . yiN,  V2N ,  2/3iv) 

(zilj  *21  >  *31, . Z1N,  Z21 V,  z3 n)-  (!) 

Figure  2  shows  the  co-ordinate  system  used  in  the  following  exposition.  ABC  is 
one  triplet  of  hydrophones  lying  on  an  equilateral  triangle,  the  plane  of  which  is 
perpendicular  to  the  tow  direction  <j>  =  0.  We  shall  ignore  the  0  dependence,  only 
considering  the  horizontal  plane  for  which  0  =  ^ .  The  desired  beam  is,  then,  steered 
towards  (6  =  <j>  =  <j>o). 
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Figure  2  Co-ordinate  system  used  in  the  report .  The  azimuthal,  horizontal  plane 
has  6  —  The  x-axis  is  defined  by  the  direction  (6  =  %,,</>  =  0). 

The  response  of  the  bearoformer,  to  a  plane  wave  coming  from  the  direction  <fi  at 
one  frequency  /,  is 

N  3,  2?t 

=  ^exp(-27r ifr3)  ^  cfy  exp(  —  fi(xkj  cos<j>  +  ykj  sin  <j>  -  cr (2) 

j-]  A.--1  c 

c  is  the  speed  of  sound,  the  r’s  and  rc’ s  are  delays  the  latter  representing  delays 
associated  with  the  individual  triplets  and  the  acs  are  the  individual  hydrophone 
weightings  ordered  on  the  triplets. 

If  we  note  that,  for  an  idealised,  non-snaking,  array,  x kj  =  xj  since  the  x  co-ordinates 
are  constant  within  a  particular  triplet: 


r(<j> ,  f)  =  Yi  exp(  27r-  ( XjCos<j>  -  cr3))  ^  akc3  exp  (—fi(ykj  sin  <j>  -  cr -*J))  (3) 

j— i  fc— i 

ajp  is  necessarily  chosen  real  (because  of  the  constraints  of  Saclant  Centre’s  present 
real  time  system). 
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It  was  found  after  some  investigation  (see  annex  A  for  proof  and  annex  E  for  limi¬ 
tations)  that  cardioid  beampatterns  are  formed  in  individual  triplets  if  a ^  =  sin  jilc 
where  filc  is  the  angular  displacement  of  the  triplet  from  the  upgoing  vertical,  and 
are  chosen  to  displace  the  hydrophone  response  to  the  origin  ’O’  in  figure  2. 

To  make  this  clearer  and  more  concrete,  we  refer  to  figure  3,  in  which  j3  is  the 
rotation  of  the  triplet 

crlj  __  r  sin^ 
cr2j  _  r  sin^3  __  7), 
cr3j  _  r  sjn^  _|_  y)^ 

where  7  is  120  degrees  (^),  and  r  is  the  radius  of  the  circumscribing  circle, 
similarly 

alJ  —  sin/3, 
a2c3  =  sin  (/?  -  7), 
cPJ  =  sin  (/3  +  7). 

It  is  found  that  the  second  sum  of  equation  3  forms  a  cardioid  shape,  at  a  given 
frequency,  which  is  almost  independent  of  the  rotation  (3.  The  magnitude  of  and 
reason  for  this  variation  are  explained  in  annexes  A  and  E.  If  we  define  this  cardioid 
response 

3  2t 

Rc(4>,f)  =  ^2<*c3  exP( — fiiVk]  sin  <p  —  CTc3)), 
k=  1  c 

which  is  approximately  the  same  for  all  cardioids,  then  equationresponse 

N  _27T  if 

r{<f>,f)  -  X!ex P( - J-{xjCos<t>  -  ct3)) 

To  steer  this  expression  we  take  cr3  =  Xkj  cos0o>  the  second  summation  can  then  be 
written  as  S(f,  </>,</>  0)  which  is  no  more  than  the  response  for  a  linear  array  steered 

tO  (j) Q. 

The  total  response  is  then 

r(f,<j>,<j>o)  =  Rc(f,(j))S{fAA  0) 

It  is  worth  noting  that  this  is  made  up  of  two  processes.  The  first  R  is  formed  on 
each  of  the  individual  triplets  (first  stage  of  processing).  This  means  each  triplet 
is  beamformed  into  a  cardioid  response.  These  ‘cardioid  hydrophones’  are  then 
themselves  beamformed  as  a  linear  array. 

-  6- 
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Figure  3  The  angle  j3  shown  for  the  array.  It  is  assumed  known  for  each  of  the 
triplets  of  the  array  allowing  the  ( x,y,z )  co-ordinates  for  the  all  hydrophones  in  the 
array  to  be  known. 


Figure  4  shows  some  2  dimensional  beampatterns,  for  9  =  f ,  for  a  range  of  azimuthal 
steer  angles.  These  plots  have  been  calculated  for  a  very  short,  10  triplet,  half¬ 
wavelength  spacing  between  triplets,  array  so  as  to  accentuate  the  effects  as  the 
look  direction  is  steered  away  from  broadside.  All  figures  have  been  normalised  to 
the  same  maximum  value  for  ease  of  comparison.  The  polar  plots  are  on  a  dB  scale 
showing  the  maximum  30dB  of  response. 

Figure  4a  shows  the  beampattern  for  the  broadside  beam.  It  is  clear  that  there  is 
only  one  significant  lobe.  This  is  of  course  the  very  result  that  we  require  (  in  the 
standard  linear  array  the  ambiguous  beam  would  be  present  pointing  towards  270 
degrees).  We  note  that  the  first  sidelobes  are  at  approximately  -13dB  (as  for  the 
ambiguous  array,  with  a  constant  shading  function) . 

Figures  4 (bed)  show  the  response  as  the  look  direction  is  steered  progressively  to¬ 
wards  forward-endfire.  As  the  look  direction  approaches  endfire  the  beampattern 
becomes  progressively  more  distorted.  It  is  clear  that  in  figure  4(d)  the  beampattern 
is  very  far  from  what  one  might  desire.  Furthermore,  the  maximum  response  of  the 
beampattern  is  no  longer  in  the  nominal  look  direction.  The  reasons  for  these  ef¬ 
fects  are  twofold.  Firstly,  as  for  the  ambiguous  array,  as  we  steer  towards  endfire  the 
mainlobe  beam  becomes  wider,  in  the  ambiguous  array  we  reach  a  point  where  the 
two  ambiguous  beams  coalesce,  in  the  directional  array,  however,  the  suppression 
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makes  the  coalesced  beam  appear  skewed.  The  second,  associated,  effect  is  more 
subtle  and  we  discuss  it  here  in  more  detail. 

Combination  of  beampatterns 

Figure  5  shows  a  diagrammatic  approach  to  understanding  how  the  Centre’s  algo¬ 
rithm  works.  In  figure  (a)  we  have  plotted  the  beampattern,  on  a  log  scale,  of  a 
single  triplet  in  blue  with  its  null  in  the  270  degree  (starboard)  direction.  On  the 
same  figure  we  have  plotted  the  response  of  a  10  element  linear,  ambiguous,  array  in 
blue.  It  has  been  steered  to  45  degrees  and  shows,  as  expected,  the  ambiguous  beam 
at  135  degrees.  The  cardioid  array  response  is  shown  in  b.  This  can  be  obtained  by 
multiplying  together  the  two  responses  of  (a).  (NB  because  the  responses  of  figure 
(a)  have  been  displayed  on  a  log  scale  some  care  must  be  taken  in  moving  from  figure 
(a)  to  (b).) 

This  understanding  of  how  the  responses  must  be  combined  is  the  first  step  towards 
explaining  the  displaced  beams  at  endfire.  In  figure  c  we  show  a  blown  up  version 
of  the  responses  for  a  beam  nominally  steered  towards  5  degrees  from  endfire,  but 
here  not  on  a  logarithmic  scale.  It  is  important  to  notice  that  the  product  of  the  two 
responses  (in  green)  does  not  have  its  maximum  at  the  maximum  of  either  of  the 
individual  responses.  This  means  that  the  overall  cardioid  array  response  maximum 
tends  to  skew  away  from  endfire  as  shown  in  figure  5(c). 

When  the  3d  algorithm  is  discussed  in  annex  B  it  will  be  seen  that  again  that 
algorithm  gives  skewed  results  towards  endfire.  There,  however,  it  will  be  prudent 
to  discuss  the  effect  in  terms  of  conflicting  constraints.  The  underlying  reason, 
however,  is  that  as  we  approach  endfire  we  progressively  lose  the  x-direction  extent, 
i.e.  the  array  appears  to  be  narrower  and  the  ability  to  resolve  to  left  and  right 
is  diminished.  At  endfire  itself  the  system  has  no  LR  extent  and  consequently  the 
system  cannot  be  expected  to  work,  regardless  of  algorithm  used. 

It  is  important  to  note  that  the  algorithm  as  expressed  above  is  a  broadband  al¬ 
gorithm  despite  the  fact  that  it  has  been  explained  and  discussed  in  terms  of  its 
individual  frequency  components.  We  can  expect,  however,  that  the  gain  achieved 
by  the  array  is  frequency  dependent.  This  is  not  surprising  since  the  directivity  in¬ 
dex  of  the  array  is  frequency  dependent.  We  shall  not  investigate  this  phenomenon 
here,  other  than  to  point  out  that  the  results  of  annex  C  can  be  used  to  quantify 
the  response  of  the  two  algorithms  detailed  in  this  report  in  an  analytic  manner. 
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Figure  4  Four  2-dimensional  beampatterns  for  various  steer  directions  in  azimuth. 
As  the  beams  are  steered  towards  endfire  the  mainlobe  becomes  progressively  distorted 
and  eventually  the  maximum  response  is  skewed  away  from  the  nominal  steer  direc¬ 
tion.  Steer  angles  are  for  (a),  (b),  (c)  and  (d)  respectively:  90  ( broadside ),  60,  SO 
and  5  degrees . 
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(a)  (b) 


Figure  5  Diagrammatic  representation  of  how  the  Centre  algorithm  works.  In 
(a)  the  ambiguous  linear  array  response  is  shown  in  blue  with  the  triplet  ‘ shading ’ 
superimposed  in  red.  (b)  shows  the  response  of  an  entire  cardioid  array  with  10 
triplets,  (c)  shows  a  blowup  of  a  10  triplet  cardioid  array  steered  to  within  5  degrees 
off  endfire.  The  red  curve  shows  the  single  triplet  response ,  the  blue  curve  shows  the 
broad  beam  formed  by  the  ambiguous  array  with  its  maximum  at  the  correct  point. 
The  product  of  these  curves  is  the  cardioid  array  response  and  is  shown  in  green. 
It  is  clearly  skewed  away  from  endfire  and  has  its  maximum  response  at  the  wrong 
point  (in  actual  fact  at  about  SO  degrees  off  endfire). 
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BACCHUS’98 


In  September  of  1998  Saclant  Undersea  Research  Centre  held  a  joint  trial  with  DERA 
(UK)  using  NRV  Alliance.  The  trial,  designated  BACCHUS’98,  had  the  primary 
objectives  of  evaluating  ambiguity-resolving  towed  arrays;  in  particular  comparing 
triplet  and  twin  array  technology.  However,  due  to  technical  problems  only  a  triplet 
array,  provided  by  DERA  was  used  to  gather  data  during  the  trial  The  array  con¬ 
sisted  of  120  hydrophones  arranged  into  40  triplets  with  a  nominal  upper  operating 
frequency  of  1780Hz. 

DERA  (UK)  also  provided  the  acoustic  source  for  the  experiment.  The  technology 
used  was  a  free-flooded  ring  (FFR).  The  typical  output  power  for  the  source  was 
220  dB  re  1/zPa.  For  the  majority  of  the  runs  carried  out  a  LFM  sweep  was  used 
between  1100  and  1800  Hz  with  a  duration  of  3.5s.  Indeed  for  all  the  data  discussed 
in  this  report  we  only  consider  waveforms  of  this  nature  and  infact  will  only  process 
a  100Hz  band  between  1450  and  1550Hz. 

The  DERA  free-flooded  ring  (FFR)  broadband  source  was  used  which  transmitted 
a  LFM  sweep  between  1100  and  1700Hz  over  3.0  seconds  followed  by  a  CW  pulse  at 
1050Hz  of  1.0  seconds  duration. 

No  detailed  information  about  conditions  or  statistical  analysis  of  data  is  given  in 
this  report.  This  will  be  deferred  to  a  subsequent  report. 

The  majority  of  runs  were  of  the  type  shown  in  figure  6  in  which  the  target  travels 
along  a  track  for  its  image  appears  to  take  it  through,  beyond  and  back  through  a 
directional  reverberation  ridge  which  had  previously  been  located  on  the  Centre’s 
SWAC  series  of  sea  trials.  For  a  conventional  array  the  target  will  become  obscured 
as  its  images  passes  through  the  ridge.  For  the  directional  array,  however,  the  ridge 
and  target  are  separated  to  port  and  starboard  so  that  the  image  remains  unobscured 
for  the  whole  run. 
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Figure  6  A  typical  run  from  the  Bacchus  Trial  of  1998.  The  main  reverberation 
ridge  running  approximately  North  South  was  used  for  many  of  the  runs  since  it  pro¬ 
vided  a  historically  well  known  source  of  directional  reverberation.  In  the  case  shown 
the  NRV  Alliance  sails  from  South  to  North.  An  echo  repeater  and  target  proceed 
on  the  v-shaped  courses  as  shown.  The  target’s  course  is  such  that  its  ambiguous 
return  appears  to  pass  into,  through  and  back  through  the  reverberation  ridge  shown. 
The  echo  repeater  is  configured  in  such  a  way  as  to  delay  the  return  transmission  so 
as  to  appear  beyond  the  reverberation  ridge.  The  Campovega  Rig  is  also  marked  for 
completeness  since  it  is  observable  in  subsequent  figures. 
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4 

Narrowband  Processing  of  Real  Data 


The  data  from  BACCHUS’98  is  inherently  broadband,  and  indeed  the  treatment 
of  section  2  is  appropriate  to  a  system  with  significant  fractional  bandwidth.  In 
this  section,  however,  we  shall  give  some  indication  of  how  the  system  works  for  a 
relatively  narrowband  data  set. 

The  original  hydrophone  data  is  filtered  to  a  band  of  100Hz  centred  at  1500Hz.  And 
then  match-filtered.  The  algorithm  of  section  2  is  used  to  beamform  the  data.  At 
the  same  time  we  take  one  hydrophone  from  each  of  the  triplets  in  the  array  to  allow 
us  to  do  traditional  beamforming  for  comparison. 

Since  the  thrust  of  this  report  is  to  detail  and  discuss  cardioid  algorithms  rather 
than  to  give  a  detailed  analysis  of  BACCHUS’98  data  we  have  investigated  only  one 
representative  run  (run  15b)  which  is  of  the  type  described  in  the  previous  section. 
Similarly  because  we  initially  wish  to  consider  and  compare  beamforming  algorithms 
no  normalization  has  been  carried  out. 

The  algorithm  of  section  2  was  used.  The  twist  angles  were  obtained  from  the  non¬ 
acoustic  data.  The  results  obtained  are  not  shown  here.  But  it  is  important  to  state 
that  the  twist  angles  were  extremely  stable  over  the  space  of  10’s  of  minutes.  This 
implies  at  least  for  the  short  array  under  investigation  here  that  the  twist  angles 
need  not  be  measured  very  often. 

When  the  beamforming  is  carried  out  it  is  over  resolved  at  1  degree  intervals.  (For 
independent  beams  the  spacing  is  about  3  degrees  at  broadside).  There  are  various 
reasons  and  consequences  attendant  on  this  decision. 

i)  because  we  are  highly  over-resolving  at  endfire  the  images  for  both  the  cardioid 
and  linear  beamforming  will  appear  to  be  ‘smeared  out’  in  the  endfire  positions  (  0, 
±80  degrees). 

ii)  we  have  obviated  the  need  to  calculate  the  3dB  crossover  points  as  we  approach 
endfire.  For  the  cardioid  array  this  is  a  non-trivial  matter  since  as  was  discussed  in 
section  2  as  we  approach  endfire  the  beams  become  skewed  and  no  longer  point  in 
the  nominal  look  direction. 

Figure  7  shows  the  result  for  one  ping  from  the  ambiguous  traditional  array.  We 
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have  plotted  both  the  port  and  starboard  responses  which  are  of  course  mirror 
images.  0  degrees  indicates  the  tow  direction  Figure  8  shows  the  same  set  of  data 
for  the  triplet  array  after  cardioid  processing.  Negative  angles  indicate  port  whilst 
starboard  angles  are  denoted  by  positive  angles. 

This  ping  shows  various  features  of  interest  which  are  worth  labouring.  The  main 
feature  is  the  curved  reverberation  ridge  which  extends  between  16  and  20km  range 
and  over  all  bearings.  This  feature  can  be  associated  with  the  long  reverberation 
ridge  shown  in  figure  6  and  discussed  at  some  length  in  the  previous  section.  After 
cardiod  processing  it  has  been  resolved  to  starboard  (in  agreement  with  figure  6!). 

At  a  range  of  approximately  13.5km  a  distinct  very  bright  bottom  feature  can  be 
discerned.  After  cardioid  processing  we  see  that  the  feature  is  found  physically  on 
the  starboard. 

At  a  range  of  24.5km  towards  forward  end-fire  a  bright  feature  can  be  seen.  This  is 
expected  to  be  the  Campovega  Rig,  and  after  cardioid  processing  this  hypothesis  is 
further  substantiated  by  the  fact  that  it  lies  to  port  (in  agreement  with  figure  6). 

One  further  feature  of  interest  is  the  other-ship  noise  which  can  be  seen  as  a  promi¬ 
nent  vertical  line  at  all  ranges.  After  processing  this  line  has  been  resolved  to 
starboard.  This  brings  out  the  interesting  point  that  because  left-right  resolution 
is  being  carried  out  on  receive  the  system  can  also  operate  in  a  passive  mode.  (As 
opposed  to  left-right  resolution  methods  which  attempt  to  resolve  on  transmit). 

Figure  9  makes  these  ideas  rather  more  concrete  and  gives  a  more  quantitative  mea¬ 
sure  of  the  left-right  suppression  achievable  with  the  beamforming  for  this  particular 
system.  It  shows  three  A-scans,  of  the  Campovega  Rig,  the  reverberation  ridge  and 
of  the  bright  bottom  feature.  The  suppression  on  the  bottom  feature  is  about  lOdB 
and  for  the  Campovega  Rig  it  is  13dB.  On  the  reverberation  ridge  a  quantitative 
measure  is  harder  to  assess  but  it  is  clear  that  the  overall  shape  of  the  ridge  has 
disappeared  (been  significantly  suppressed)  on  the  port  side  of  the  display. 
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Figure  7  B-scan  for  non- directional  array.  The  display  has  been  plotted,  over  the 
whole  range  of  angles ,  both  port  and  starboard  to  facilitate  comparison  with  figure  8. 
The  boxed  bright  spots  at  a  range  of  approximately  25km  is  the  Campovega  Rig.  The 
boxes  at  13.5km  highlight  a  bright  bottom  feature.  The  vertical  line  at  zb  180  degrees 
is  other  passive  ship  noise . 
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Figure  8  B-scan  for  directional  array.  In  this  display  it  is  clear  that  the  Campovega 
Rig  is  on  the  port  side  and  the  reverberation  ridge  and  the  other-ship  noise  have  been 
resolved  to  starboard  in  agreement  with  figure  6. 
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Bearing  in  Degrees 


Figure  9  A -scans  of  three  representative  points  taken  from  the  B-scan  of  figure  8. 
The  upper  trace  exhibits  the  Campovega  Rig  which  shows  a  suppression  of  approxi¬ 
mately  12dB  between  port  and  starboard.  In  the  middle  scan  the  reverberation  ridge 
can  be  seen  to  have  been  suppressed  on  the  port  side.  In  the  bottom  scan  the  bright 
bottom  feature  can  be  seen  to  have  been  suppressed  by  approximately  lOdB.  (Note 
the  ordinate  scale  is  different  for  the  three  displays) 
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5 

Conclusions 


The  main  purpose  of  this  report  is  the  documentation  of  the  proposed  cardioid  algo¬ 
rithm  which  is  to  be  used  by  the  Saclant  Undersea  Research  Centre.  Consequently, 
some  effort  has  been  expended  in  detailing  the  algorithm  and  quantifying  some  of 
its  particular  properties.  The  appendices  have  been  used  to  tabulate  the  analytic 
expressions  for  SNR  and  left-right  suppression. 

We  have  taken  the  opportunity,  however,  to  describe  and  compare  a  ‘3d-algorithm’. 
We  have  shown  how  the  ‘3d-algorithm’  gives  superior  left-right  suppression  but  the 
‘Centre  algorithm’  gives  superior  SNR  performance  in  uncorrlated  noise. 

It  is  recommended  that  this  work  should  be  further  verified  and  extended  by  ap¬ 
plication  to  more  data.  In  particular  the  application  of  the  cardioid  algorithms 
to 

1.  investigate  the  difference  in  left-right  reverberation  statistics.  A  directional 
array  gives  an  investigator  a  unique  opportunity  to  compare,  for  instance,  the 
difference  between  the  characteristics  of  upslope  and  downslope  reverberation 
statistics. 

2.  Being  broadband  the  algorithm  allows  the  investigation  of  not  only  left-right 
comparisons  but  also  the  processing  of  different  bands.  For  instance  the  Bac¬ 
chus  data  (section  3)  although  only  600Hz  wide  is  known  to  show  very  different 
performance  for  the  upper  and  lower  bands  [4].  It  is  highly  recommended  that 
these  effects  should  be  further  investigated  and  quantified  as  an  initial  indica¬ 
tor  for  the  Mercury  trial  of  1999. 
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Annex  A 

Real  Weights 


In  section  2  we  used  a  real  weighting  on  each  of  the  individual  triplet  hydrophones 
which  was  proportional  to  the  projection  of  the  hydrophones  position  vector  onto  the 
horizontal.  We  did  this  because  of  limitations  in  the  Centre’s  real  time  beamformer. 
In  this  annex  we  discuss  the  approximations  inherent  in  this  approach. 

If  we  consider  a  single  triplet  and  impose  the  following  constraints 

wHci  =  1  (4) 


where  ci  is  the  steering  vector  defining  the  look  direction  and  w  is  the  weight  vector 
that  we  wish  to  calculate.  For  the  triplet  this  lies  in  the  horizontal,  in  the  plane 
of  the  triplet  and  the  actual  direction  defines  whether  our  beamformed  cardioid  is 
looking  left  or  right.  We  also  impose 

whc2  =  0  (5) 

C2  is  the  steering  vector  for  the  opposite  direction  to  the  look  direction  and  as  such 
defines  the  null  in  our  beampattern. 

These  two  constraints  can  be  written  more  succinctly  as 

CHw  =  b  (6) 


with 

C  =  [C1C2]  and  b  —  [1,0]T 


Furthermore,  let  us  define, 


WHW  =  1,  CiHd  =  1.0,  C2HC2  =  1.0, 


for  mathematical  convenience. 

Assuming  that  noise  is  uncorrelated  between  the  triplet  hydrophones,  the  optimum 
set  of  weights,  i.e.  which  maximises  the  SNR,  is  given  by. [6] 

wopt  =  CfCnC)"1!*.  (7) 
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Which  after  some  rearrangement  becomes: 

wopt  oc  ci  -  a* c2,  (8) 

where  a  —  Cj*C2  is  the  correlation  between  the  vectors  Ci  and  C2.  It  is  interesting  to 
note  that  if  Ci  and  C2  are  orthogonal  the  optimum  weight  is  no  more  than  the  steering 
vector  (the  usual  case  in  beamforming).  We  shall  return  to  these  considerations  again 
in  annex  C  where  this  term  plays  a  crucial  role  in  the  signal  to  noise  ratio  gain  of 
the  system. 

At  this  stage  it  is  helpful  to  consider  the  individual  terms  of  the  above  expression 


wopt  = 


27razsin/9 

e27r ai  sin  j3 

1 

e~2nai  sin((3 —7) 

a* 

e27raisin(/3— 7) 

V3 

e—27r  ai  sin(/?+7) 

~  A 

e27rai  sin(^-f  7) 

For  the  triplet  we  can  write 


a*  —  l/-47rcusin/?  g— 47tcw  sin(/3— 7)  _|_  e-47rmsin(£-7)^ 

a  - 


(9) 


If  the  parameter  a,  which  is  the  radius  of  the  circumscribing  radius  of  the  triplet  in 
units  of  wavelength,  is  small  we  can  make  expansions  to  first  order  in  a. 

a*  —  -(3  —  47ra*(sin/3  +  sin(/3  —  7)  +  sin(/3  —  7)).  (10) 

3 


Now,  we  note  that,  from  annex  F 

sin/3  +  sin(/3  -  7)  +  sin(/3  +  7)  =  0.0 


so  that 

then 

and  finally 


=  1.0 


Wopt  = 


V3 


e2Jroisin/?(1  _  (i  +  4tt<m  sin  /?)) 

_  (1  +  47rai sin (J3  -  7)) 
e27T<nsm(/3+7)(l  -  (1+  47raisin(/3  +  7)) 


wopt  oc 


sin  j3 

sin{p  -  7) 

© 

sin(/3  +  7) 

-2tt ai  sin (3 


e~2irai  sin  (j3 — p) 
g-27 xai  sin (£+7) 


where  ©  is  the  element-by-eiement  product  operator. 


(ii) 


(12) 
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We  have,  then,  that  the  optimum  weights  for  the  triplets  with  the  given  constraints 
have  the  effect  of  delaying  the  signal  to  the  centre  of  the  circumscribing  circle  of 
the  triplet  (given  by  the  exponential  phase  terms)  and  real  weights  which  are  in  the 
following  ratio 

K1  _  lw2|  _  lw3|  Q3) 

sin  (3  sin(/3  —  7)  sin(/3  +  7) 

At  this  point  it  is  worth  stressing  that  although  the  above  calculation  has  been 
carried  out  as  a  narrowband  calculation  the  final  result  for  the  optimum  weights 
for  a  broadband  signal  is  made  up  of  a  set  of  delays  and  real  weights  which  are 
independent  of  frequency,  assuming  that  the  small  a  expansion  is  valid  for  at  all 
points  in  the  band.  Consequently  the  result  holds  across  any  band  of  interest.  The 
observant  reader,  however,  will  have  noticed  that  the  small  term  expansions  carried 
out  above  have  not  been  carried  out  completely  i.e.  we  have  retained  all  powers 
of  a  in  the  delay  term  of  equation  11.  The  reason  for  this  is  that  we  wish  to  find 
an  explanation  and  justification  for  the  expression  for  the  triplet  weights  given  in 
section  2. 
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Annex  B 

Alternative  Narrowband  Approach 


In  the  main  body  of  the  text  we  have  constrained  our  prospective  algorithm  to  be 
composed  of  two  stages,  the  first  an  initial  beamforming  at  the  triplet  level.  The 
second  a  Traditional’  steering  of  the  array  using  these  new  ‘directional  hydrophones’. 
This  approach  is,  to  a  certain  extent,  imposed  on  us  by  the  structure  of  Saclant 
Centre’s  real-time  system. 

If  we,  however,  consider  the  entire  array  as  a  3d  beamformer  and  limit  ourselves  to 
narrowband  systems  the  choice  of  algorithms  becomes  substantially  widened. 

In  this  section  we  shall  simply  state  an  obvious  narrowband  algorithm.  We  shall, 
however,  return  to  it  in  more  detail  when  we  are  considering  some  of  the  signal-to- 
noise  ratio  properties  of  the  two  algorithms. 

With  the  freedom  offered  by  the  full  three  dimensional  nature  of  the  array  the 
most  obvious  approach  is  to  impose  the  following  constraints,  using  the  notation  of 
annex  A: 

cfw  =  1 
c^w  =  0 

with,  ci,C2  representing  the  vectors  which  point  in  the  look  direction  and  the  null 
direction  respectively.  Here  we  are  interested  in  the  situation  where  the  null  direction 
is  in  the  ambiguous  direction  to  the  look  direction  and  in  the  same  plane.  As  usual 
we  define  the  constraint  vectors  such  that  |ci|2  =  1  and  |c2 1 2  =  1*  Now,  however, 
these  vectors  contain  3 N  elements  where  N  is  the  number  of  cardioids.  We  combine 
the  two  vector  constraints  into  one  matrix  constraint: 

CH  =  b  (14) 

where 

C  =  [Clca],  b  =  [10]H 

If  we  take  the  pseudo-inverse  we  find,  after  some  arrangement,  the  following  solution 
for  w  analagous  to  that  given  in  annex  A: 

W=- — T^(C1  -  a*c2)  (15) 

1  —  |  a\z 

where  a  =  c^C2  is  the  correlation  between  the  vectors  Ci  and  c2. 
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Annex  C 

Signal  to  Noise  Ratio  Considerations 


In  this  annex  we  will  point  out  some  of  the  signal  to  noise  properties  of  the  algorithms 
described  in  the  body  of  the  main  text  and  the  3-dimensional  approach  of  annex  B. 
For  both  of  the  algorithms  we  shall  consider  the  signal  gain,  the  uncorrelated  noise 
gain  and  the  correlated  noise  gain. 

We  begin  by  considering  the  performance  of  the  two  algorithms  in  uncorrelated  noise. 
This  is  the  noise  that  we  associate  typically  with  electronic  noise  in  the  system. 
There  is  also  a  component  of  flow  noise  which  can  be  modelled  as  uncorrelated  even 
within  a  triplet. [7] 

Next  the  performance  of  the  two  algorithms  is  analysed  in  correlated  noise  which 
provides  a  model  of  sea-noise.  It  is  worth  noting  that  the  performance  of  the  algo¬ 
rithm  in  ‘directional  noise’  is  best  considered  in  terms  of  left-right  suppression  which 
is  addressed  in  annex  D. 


Signal  to  Uncorrelated  Noise  for  Centre  Algorithm 

For  the  Centre  algorithm  we  shall  assume  that  the  triplets  are  not  rotated.  It  is 
shown  in  annex  E  that  the  overall  response  of  the  triplet  remains  fairly  constant 
with  rotation  so  that  the  results  found  here  are  general.  We  shall  only  consider 
signals  arriving  in  the  horizontal  plane.  This  is  of  course  an  approximation  since 
in  a  real  sonar  system  we  may  expect  returns  from  above  or  below  horizontal.  We 
shall  consider  the  signal  and  the  noise  output  powers  separately. 


Signal  Component 


For  a  general  plane-wave  signal  on  the  horizontal  coming  from  an  azimuthal  angle 
of  (j)  the  response  is  given  by: 


1 

7K 


Nc  3 

Ec,E  Cf  ex.p(2ni(asinryi  —  Xj  cos  <j>  —  yi  sin  <fi)) 
j=l 


(16) 


where  Cj  represents  the  complex  weight  for  the  steering  in  azimuth,  Cf  is  the  real 
cardioid  weight,  7 i  is  the  rotation  of  the  particular  hydrophone  i,  within  the  triplet, 
from  the  normal.  It  is  clear  after  a  little  thought  that 

yi  =  a  sin  7;. 
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And  so  after  some  manipulation  we  can  write  the  response  as 
_ 3 

\fN~c  y:  sin  7 i  exp(27T ai  sin  7^(1  -  sin  </>)). 
i 

Since  we  assume  that  no  rotation  has  taken  place  sin  71  =  0,  sin  72 
sin 73  =  — So  that  we  can  write  the  output  as 

i \J 3 Nc  sin  a/3 ttu  ( 1  —  sm(j>) 


(17) 


^  and 


and  the  output  signal  power  is 

Output  Signal  Power  =  3NC  sm2(V3wa(l  -  sin  </>)).  (18) 


With  broadside  at  </>  = 

Uncorrelated  Noise  component  for  Centre  Algorithm 


Calculation  of  the  noise  output  power  is  similarly  straightforward  although  it  can 
only  be  expressed  in  a  statistical  manner.  The  expectation  of  the  output  power  is 


(19) 


where  rtij  represents  a  noise  sample  on  the  ith  hydrophone  of  the  jth  triplet. 

=  E  E  E  E 

i  V  j  j' 

=EEi^mr- 

*  3 

But  on  using  the  explicit  values  for  these  parameters  and  using  the  results  of  ap¬ 
pendix  F  we  obtain 

Noise  Output  Power  =  1.5, 


and  finally  we  can  write 

Output  Signal  to  Noise  Ratio  =  2 Nc  sin2(\/37ra(l  -  sin  <j>))  (20) 


This  result,  as  it  stands  is  far  from  self-explanatory  so  some  comments  are  worth- 
while. 

(i)  This  result  has  been  calculated  for  two  hydrophones  on  the  horizontal.  When 
their  horizontal  spacing  is  ^  the  Output  Signal  to  Noise  Ratio  (which  in  this  case  is 
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the  gain  of  the  system)  for  a  signal  from  broadside  is  2 Nc.  This  compares  with  the 
3 Nc  we  would  usually  expect  for  a  system  with  3 Nc  hydrophones.  This  is  because, 
in  effect  one  of  the  hydrophones  has  been  turned  off. 

(ii)  as  a  decreases  the  signal  to  noise  ratio  also  falls.  The  reason  for  this  is  due  to 
the  decrease  in  signal  rather  than  any  change  in  noise  output  power.  As  an  order  of 
magnitude  indication  if  we  assume  that  a  =  0.05 A  as  a  typical  value  then  the  output 
SNR  at  broadside  is  2 Nc  0.27.  A  loss  of  almost  6dB  over  the  ^  case. 

(iii)  As  the  look-direction  is  steered  from  broadside  the  output  signal  to  noise  ratio 
decreases.  If  we  take  the  loss  of  output  signal  to  noise  ratio  at  endfire  compared  to 
that  at  broadside  we  get; 


SNRioss  = 


sin2  (y/Sna) 
sin2(2\/3 7ra) 


(21) 


For  small  values  of  a  this  ratio  becomes  a  factor  of  a  quarter  or  in  otherwords  there 
is  a  loss  of  6dB  as  we  steer  from  broadside  to  endfire. 


(iv)  It  is  worth  noting  that  the  expression  20  above  when  plotted  out  actually  has 
the  cardioid  shape  over  the  entire  range  of  </>.  In  practice,  of  course,  only  one  half  of 
the  range  of  angles  is  used  since  an  alternative  set  of  weights  would  be  used  to  look 
on  the  opposite  side  of  the  array. 

Signal  to  Uncorrelated  Noise  for  3d  Algorithm 

We  saw  in  annex  B  how  the  optimum  weight  vector  for  the  3d  algorithm  can  be 
written  as 

wocci  -  a*C2 


We  ignore  the  constant  of  proportionality  in  what  follows  since  we  are  ultimately 
interested  in  ratios. 

w  is  the  weight  vector  matched  to  a  signal  coming  from  the  direction  given  by  Ci. 
A  signal  of  unit  power  can  be  represented,  therefore,  as  y/Mci  where  M  =  3 Nc  is 
the  total  number  of  sensors.  The  output  signal  power  then  is  given  by 

P  =  M|w?ci|2  =  M(  1  -  |a|2)2. 

For  such  a  system  the  output  uncorrelated  noise  power,  covariance  matrix  Rjy  =  / 
is 

whRnw  =  wHIw  —  wHw  =  1  —  |ai|2.  (22) 
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Then  the  output  signal-to-noise  ratio  is 

SNRout  =  M(  1  -  |of |2) 


(23) 


It  is  interesting,  then,  to  ask  the  questions  how  does  this  algorithm  perform  at 
broadside  and  as  we  steer  towards  endfire.  We  wish  to  do  this  ultimately  of  course 
with  reference  to  the  Centre  algorithm  discussed  above. 

At  broadside  then  we  can  write,  expanding  up  to  powers  of  a2 

M(1  —  |a|2)  =  3iVc(87r2a2) 

where  we  have  used  the  results  of  annex  F. 

But,  referring  to  equation  23  and  expanding  it  to  the  second  order  in  a  we  see  that 
the  signal  to  noise  ratio  at  broadside  is  equal  for  the  two  algorithms.  This  is  as 
expected  since,  at  broadside,  the  3d-algorithm  forms  a  cardioid  on  each  triplet. 

Away  from  broadside  we  can  write  generally  that 

3 a  —  c^C2  =  exp(4:7T  sin  71  sin  <f>  +  47r  sin  72  sin  <j>  +  47r  sin  73  sin  </>)  (24) 

Again  using  expansions  up  to  second  order 

a  «  1  —  47 r2a2  sin2  <fr 

so  that,  the  uncorrelated  SNR  is 

M(  1  -  |a|2)  «  3Nc(87r2 a2)  sin2  4 > 

It  is  clear  then  that  the  3d-algorithm  has  an  inferior  performance  with  regard  to 
uncorrelated  signal-to-noise  ratio  at  all  angles  steered  off  broadside. 

Performance  of  Center  Algorithm  in  Correlated  Noise. 

When  we  come  to  consider  how  algorithms  behave  in  the  presence  of  correlated 
noise  expressions  rapidly  become  complicated  and  unilluminating.  However,  we  shall 
consider  the  Centre  algorithm  in  some  detail  since  it  represents  a  crucial  measure 
for  real  experimental  systems. 

We  assume  the  narrowband  approximation  for  the  covariance  matrix: 

/1  /  f\ 

Rn  =  /  1  /  (25) 

\f  f  lJ 


SACLANTCEN  SR-329 


-  26 


SACLANTCEN  SR-329 


where 


/ 


sin(27ra) 

2iva 


and  a  =  y/Z a  is  the  distance  between  hydrophones.  The  noise  output  power  then 
for  one  triplet  is  wH R^w  which  can  be  written  as 


1  1  1  U. 
1  1  1 


(1  —  f)wH  Iw  +  fwH 
The  second  term  can  be  shown  to  be  2a2f.  Consequently, 


/  111 

V 


(26) 


(1 


2) 


Signal  to  correlated  noise  power  per  triplet  =  ■  — 17  9 rj.  (27) 

1.5(1  — /)  + 

But,  taking  series  expansions  to  a2, 

sin(27ra)  ^  3(27ra)2 

27ra  3! 


and 


(i  -  /)  = 


(27ra)2 


Now  if  we  ignore  the  intertriplet  correlation  then  we  can  just  scale  this  result  by  the 
number  of  triplets  Nc. 

3  n 

Signal  to  correlated  noise  power  for  array  =  -~^(1  —  sirup)  (28) 

for  small  a. 

In  reality  of  course  the  overall  signal  to  noise  ratio  at  the  output  is  an  appropriately 
weighted  combination  of  correlated  and  uncorrelated  SNR’s. 

Performance  of  3d-Algorithm  in  Correlated  Noise 

The  analysis  of  the  correlated  noise  for  the  3d-algorithm  is  very  similar  to  that  for 
the  Centre  algorithm  and  it  is  easily  shown  that  after  only  including  terms  up  to  a2 
that  the 

Signal  power  =  3Nc(Sn2a2  sin2  <p).  (29) 


that 


Correlated  Noise  power  =  87r2a2  sin2  <p. 


Signal  to  correlated  noise  power  for  array  =  3 iVc. 


(30) 

(31) 
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This  final  expression  needs  some  comment.  It  appears  that  the  array  is  now  giving  a 
gain  of  37VC  i.e.  the  gain  that  would  have  been  achieved,  against  uncorrelated  noise, 
if  we  had  used  all  of  the  hydrophones  in  a  linear  configuration  but  now  we  have 
gained  the  ability  to  discriminate  between  left  and  right.  We  have  seen  however,  that 
both  the  algorithms  discussed  have  relatively  poor  performance,  against  uncorrelated 
noise.  The  choice  of  whether  to  use  one  of  the  algorithms  discussed  here  or  to  use  a 
linear  array  subset  (which  is  always  available  by  taking  one  element  from  each  triplet) 
is  always  open  to  the  sonar  operator.  However,  the  actual  operating  conditions  under 
which  a  triplet  algorithm  should  or  should  not  be  used  is  still  to  be  investigated. 

Also  in  any  comparison  between  directional  and  Ordinary1  processing  it  is  important 
to  realise  that  the  two  approaches  differ  at  a  very  basic  level.  Simplistically,  for 
ordinary  processing  we  add  signals,  for  the  directional  processing  discussed  above 
signals  are  subtracted.  In  the  preceding  discussions  this  has  caused  no  problems 
because  the  approach  has  been  analytic.  However,  in  reality,  such  systems  will  be 
limited  by  the  number  of  bits  being  used  in  the  calculations. 
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Annex  D 

Left-Right  suppression 


We  can  easily  calculate  the  expected  left-right  suppression  for  the  two  algorithms 
which  have  previously  been  called  the  c3d-algorithm’  and  the  ‘Centre-algorithm’. 
When  we  talk  about  left-right  suppression  here  we  are  interested  in  suppression  of 
a  signal  assuming  that  it  can  be  modelled  by  a  plane  wave. 


Centre  Algorithm 

We  can  use  equation  18  in  annex  C  to  calculate  the  left-right  suppression  for  the 
Centre  algorithm. 


The  output  power  for  a  signal  with  a  look  direction  given  by  an  azimuthal  angle  of 
4>  is 

3NC  sin2(\/37ra(l  —  sin  </>)), 

the  output  power  for  a  signal  from  the  ambiguous  direction  (180  +  </>)  is 

37VC  sin2  ( VSna  ( 1  +  sin  </>)) 


So  that  the  left-right  suppression  is  approximately, 


(1  —  sin  (f>)2 
(1  +  sin  (j>)2 


7T 

2* 


(32) 


The  dashed  lines  of  figure  10  shows  the  theoretical  left  right  suppression  as  a  function 
of  (j).  It  clearly  has  the  correct  overall  form  with  suppression  falling  from,  nominally 
an  infinite  value  at  broadside  to  1  at  endfire  where  no  left-right  discrimination  is 
possible  for  the  cardioid. 


The  crosses  show  the  values  taken  from  experimental  data  taken  from  trial  BAC- 
CHUS’98.  40  pings  have  been  analysed  choosing  objects  which  appear  point-like  on 
a  b-scan  display.  The  suppression  between  right  and  left  is  measured  and  displayed 
as  a  star  on  figure  10.  A  dense  cluster  of  returns  can  be  seen  at  140  degrees.  This 
is  the  Campovega  Rig,  which  for  the  40  pings  investigated  is  always  visible  with  a 
large  signal  to  noise  ratio,  It  is  important  to  note  that  the  dashed  line  of  figure  10 
represents  a  maximum  suppression.  Reasons  why  this  may  not  be  achieved  include; 


(i)  physical  mismatches  in  the  equipment  such  as  errors  in  the  roll  sensors,  phase 
and  amplitude  errors  in  the  hydrophones,  bandwidth  effects  (we  treat  100  Hz  band 
as  narrow  band). 
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(ii)  Channel  effects.  Any  physical  mechanism  which  stops  the  return  from  a  point  like 
object  looking  like  a  plane  wave  at  the  array  may  degrade  the  performance  (although 
it  is  conceivable  that  they  will  improve  the  performance  in  some  situations). 

(iii)  If  another  object  really  can  be  found  in  the  ambiguous  direction  the  suppression 
may  be  significantly  less  than  expected.  By  the  same  argument  the  maximum  sup¬ 
pression  that  can  be  expected  even  in  a  noise-limited  environment  is  no  more  than 
the  signal  to  noise  ratio  since  even  when  the  ambiguous  return  is  pushed  below  the 
noise  the  ambiguous  return  will  be  equivalent  to  the  noise  level. 

3d- Algorithm 

Theoretically  the  left  right  suppression  is  always  infinite,  since  the  constraint  in  the 
ambiguous  direction  is  chosen  to  be  zero.  In  reality  of  course  there  will  be  leakage 
due  to  errors  (see  annex  E)  and  the  signal  will  not  come  entirely  from  the  horizontal, 
or  indeed  from  one  particular  azimuthal  direction. 

In  figure  11  the  measured  suppression  is  shown  for  the  3d-algorithm.  The  notation 
used  is  as  for  the  Centre  algorithm.  Here  we  have  retained  the  theoretical  dashed 
curve  for  the  Centre  algorithm  to  aid  the  reader.  It  is  clear  that  the  measured  results 
often  exceed  the  Centre  algorithm  indicating  that  the  3d-algorithm  exhibits  superior 
performance.  Again  the  Campovega  Rig  returns  can  be  seen  forming  a  cluster  at 
140  degrees  but  now  showing  a  suppression  of  about  20dB. 

Figure  12  is  included  for  completeness,  and  should  ideally  be  compared  with  figure  9. 
The  suppression  of  the  bottom  feature  and  the  Campovega  Rig  is  such  that  no 
indication  of  the  strong  features  remain  above  the  noise  floor. 
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Bearing 


Figure  10  Theoretical  and  Experimental  values  for  Left- Right  suppression  for  the 
Center  algorithm  is  used.  Only  small  targets  with  large  SNR ’s  were  investigated  over 
40  pings  of  data.  The  crosses  show  the  experimentally  measured  values  of  SNR.  The 
dashed  line  indicates  the  theoretical  maximum 1  given  by  equation  32 
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Figure  11  Experimental  values  for  Left- Right  suppression  for  the  3d  algorithm. 
The  dashed  line  shows  the  theoretical  limits  for  the  Centre  algorithm  it  is  clear  that 
for  a  range  of  look  directions  the  3d  algorithm  is  superior  in  terms  of  suppression. 
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Bearing  in  Degrees 

Figure  12  A  scans  of  three  representative  points  for  the  same  data  set  as  used  in 
figure  9.  The  upper  display  shows  the  results  for  the  Campovega  Rig ,  the  middle  scan 
shows  the  Reverberation  ridge  and  the  bottom  scan  shows  a  bright  bottom  feature. 
In  distinction  to  figure  9  the  suppression  of  the  objects  is  down  into  the  noise  floor 
regardless  of  angle ,  This  is  the  same  effect  as  has  been  exhibited  more  quantitatively 
in  figure  11.  (Note  the  change  of  ordinate  scale  between  the  three  displays). 
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Annex  E 

Error  Considerations 


The  first  ‘error’  that  we  shall  address  is  that  of  twist  oil  the  sensors  on  the  ar¬ 
ray.  It  was  shown  in  annex  A  that  the  algorithm  described  in  section  2  was  an 
approximation.  The  consequence  of  the  approximation  is  that  for  any  particular 
triplet  the  response,  although  always  being  a  maximum  in  the  horizontal,  will  vary 
in  magnitude  with  twist  angle  even  if  this  value  is  known  exactly.  If  this  variation 
were  large  the  effect  on  the  beamforming  could  be  serious  since  each  triplet  would 
have  a  different  gain.  We  test  this  effect  ‘empirically’  by  calculating  the  beampat- 
terns  for  a  single  triplet  for  a  range  of  twist  angles  for  a  system  operating  at  3400 
Hz  and  a  triplet  with  radius  0.027m.  Figure  13  shows  the  response  of  a  triplet  as 
the  twist  angle  (/?  in  figure  3)  is  varied.  The  response  is  taken  on  the  horizontal 
and  is  displayed  in  dB.  Only  angles  between  0  and  60  degrees  are  shown  since  the 
response  repeats  every  60  degrees.  It  is  clear  that  the  response  variation  is  small 
(about  0.4dB).  For  lower  frequencies  the  variation  is  expected  to  be  less  (a  would 
be  smaller  in  equation  11). 

However,  in  the  real  system,  we  also  expect  an  error  to  be  associated  with  the  twist 
angle  (we  use  fi  but  the  actual  value  is  /?').  The  reasons  for  this  are  clear,  there  may 
be  an  error  in  the  twist  sensor  itself  or  since  we  interpolate  between  twist  sensors  to 
find  the  twist  at  individual  triplets  an  extra  error  may  be  introduced. 

We  ask  the  question  then,  what  is  the  response  in  the  horizontal  in  the  null  and  the 
look  direction  for  an  individual  triplet  rotated  to  /3f  if  the  actual  rotation  is  /3.  It  is 
important  to  note  that  these  are  not  necessarily  the  minimum  or  maximum  response 
for  the  whole  cardioid  response. 


The  actual  steering  vector  for  a  signal  coming  from  the  nominal  null  direction  is, 
ignoring  normalization  terms, 


_  ^2irai  sin  /?  ^2-n  ai  sin(/?— 7)  ^2tv ai  sin(/3+7)  j 


(33) 


The  weights  used  are 

w  =  [sin  /3'e2”aismP'  -  7)e2™;5in(/?'-7);  sin (p'  (34) 

Consequently,  the  response  is,  after  taking  note  of  the  fact  that  a  is  small  and  making 
the  appropriate  expansions 
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Twist  on  Cardioid 

Figure  13  The  maximum  response  of  a  triplet  after  using  the  algorithm  of  section  2 
plotted  against  known  twist  angle  (in  degrees).  The  maximum  response  is  always  in 
the  horizontal.  The  operating  frequency  is  taken  at  3400Hz  (as  opposed  to  the  1500Hz 
of  section  4)  for  a  triplet  radius  of  0.027m . 


wHs  oc  sin  (sin /3  —  sin/?')  T  sin^'  —  7)(sin(/3  —  7)  -  sin(/3/  —  7)) 

+  sin {/3'  -I  7) (sin(/3  +  7)  -  sm(/3'  +  7))  (35) 

and  using  the  results  of  section  F 

wHs  a  (1 . cos(/3  --  /?'))  (36) 

and  if  the  error  is  small 

|wHs|2  oc  (/?  -  /3')4  (37) 

In  the  nominal  look  direction  a  similar  treatment  gives 

wHs  oc  (1  +  cos(^  -  fdf))  (38) 

We  make  no  further  comment  on  these  analytic  expressions  at  the  moment >  other 
than  to  note  that  empirically  the  errors  associated  with  the  twist  sensors  seem  to 
be  relatively  unimportant 

Phase  errors 
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We  now  carry  out  a  similar  treatment  for  phase  errors.  These  errors  may  have 
two  cause:  system  errors  (e.g.  badly  matched  hydrophones)  or  signal  processing 
effects.  In  the  broadband  processing  the  phases  are  applied  as  delays  produced  by 
an  interpolation  filter,  discretization  errors  in  the  delays  will  equate  to  phase  errors 
in  the  frequency  domain.  We  expect  the  second,  signal  processing,  effect  to  be 
dominant  so  we  shall  associate  the  errors  with  the  weights  applied. 


We  assume  phase  errors  (5i,^2?^3)  on  the  hydrophones  of  the  triplet.  These  vary 
from  sample  to  sample  and  are  zero  mean.  We  also  expect  that  each  is  uniformly 
distributed;  although  this  will  not  be  used  in  what  follows.  The  steering  vector  is 


s 


r  27rmsin/3  2nai  sin(/3— 7) 
[e  5 e  5 


2-nai 


£sin(/3+7)j 


(39) 


the  weights  applied  are 

w  =  [sin£e2™isin^+l5l,sin(/3  -  7)e2’rai(sin(/3“7)fc,sin(^  +  7)e2^sm(/3+7)+^]  (40) 

when  a  is  small 

wHs  oc  sin  f3Si  -f  sin(/3  —  7)^2  +  sin(/3  +  7)^3  (41) 

but  the  Ss  are  random  variables  with  zero  mean  so  we  should  deal  with  expectations 

£(whs)  -  0  (42) 

but  more  importantly,  if  we  assume 

E(8iSj)  =  a28ij  (43) 

where  Sij  is  the  kronecker  delta.  After  using  results  from  annex  F,  we  find 

U(|wHs|2)  =  1.5<r2  (44) 
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Annex  F 

Some  Useful  Expressions 


For  ease  of  use  we  state  some  of  the  formulae  which  have  been  used  at  various  points 
throughout  this  report.  We  are  concerned  with  expressions  in  which  three  terms  are 
summed,  one  from  each  of  the  sensors  on  the  triplet.  We  shall  use  the  same  notation 
as  throughout  the  report  where  7  is  120  degrees  (  ^  radian  ),  the  angular  separation 
between  the  sensors  in  the  triplet. 

Some  tedious  algebra  gives  us 

Result  1 


sin  (3  sin  ft  +  sin(/3  —  7)  sin^  —  7)  +  sm(p  +  7)  sin {ft  +  7)  =  -cos(ft  —  ft)  (45) 


If  we  set  ft  =  j3  we  get 

Result  2 

sin2  P  +  sin2(/3  —  7)  +  sin  2(/3  +  7)  =  -  (46) 

Next  consider  the  expression 

ei(3  +  ei{P- 7)  +  e*(/?+7) 


we  can  rearrange  this  as 

eij3{l  +  e~^  +  e^)  =  e^(l  +  2cos7) 

but  cos  7  =  —  |  so  that  we  can  write,  after  taking  the  real  and  imaginary  parts: 

Result  3 


Result  4 


sin  P  +  sin(/3  —  7)  +  sin(/3  +  7)  =  0.0 


(47) 


cos  p  +  cos(/3  —  7)  +  cos  (p  +  7)  =  0.0 


(48) 
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Abstract 


It  is  well  known  that  conventional  active  towed  arrays  suffer  from  ‘The  Ambiguity  Problem’  in  which  it  is 
impossible  to  distinguish  returns  from  port  or  starboard.  Although  operational  methods  for  overcoming 
this  problem  exist  they  are  time-consuming  and  cannot  be  used  for  a  single  ping.  One  proposed  method  of 
coping  with  these  difficulties  is  to  use  a  triplet  array  in  which  the  direction  of  arrival  of  a  signal  can  be 
ascertained  on  a  single  ping.  This  report  details  a  broadband  algorithm  which  has  been  developed  and  is 
currently  being  used  at  the  Saclant  Undersea  Research  Centre.  The  algorithm  was  specifically  designed  to 
be  used  on  a  delay-type  beamformer  in  which  shading  cannot  be  applied  at  the  beamlevel.  4  second 
algorithm  employing  the  full-flexibility  of  a  3d  beamformer  is  also  described  and  the  two  algorithms  are 
compared  for  their  performance  in  terms  of  left-right  suppression  and  SNR  for  which  analytical 
expressions  are  computed.  It  is  shown  that  both  of  the  algorithms  have  desirable  but  different  features. 
The  algorithms  are  applied  to  data  taken  from  the  trial  BACCHUS’98  and  the  results  obtained  are  used  to 
help  validate  the  theoretical  results  and  to  give  an  indication  of  the  performance  which  is  achievable  for 
the  triplet  system.  We  show  that  for  the  Centre  algorithm  a  typical  achievable  left-right  suppression  is 
greater  than  15dB. 
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